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A Bose-Einstein condensate exhibiting a nontrivial phase induces an artificial magnetic field in immersed 
impurity atoms trapped in a stationary, ring-shaped optical lattice. We present an effective Hamiltonian for 
the impurities for two condensate setups: the condensate in a rotating ring and in an excited rotational state in 
a stationary ring. We use Bogoliubov theory to derive analytical formulas for the induced artificial magnetic 
field and the hopping amplitude in the limit of low condensate temperature where the impurity dynamics is 
coherent. As methods for observing the artificial magnetic field we discuss time of flight imaging and mass 
current measurements. Moreover, we compare the analytical results of the effective model to numerical results 
of a corresponding two-species Bose-Hubbard model. We also study numerically the clustering properties of 
the impurities and the quantum chaotic behavior of the two-species Bose-Hubbard model. 
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I. INTRODUCTION 

A prominent application of ultracold atoms trapped in opti- 
cal lattices is that of a quantum simulator of condensed mat- 
ter models. The advantage of such a simulator is that the 
Hamiltonian defining the condensed matter model can be im- 
plemented almost perfectly, and the relevant parameters, such 
as interaction strength or lattice geometry, can be tuned in 
a wide parameter range [1,2]. In contrast, for a given con- 
densed matter system the Hamiltonian governing its dynam- 
ics is often difficult to find. Examples of effects which can be 
simulated with neutral atoms in optical lattice are the quantum 
Hall effect in two dimensions [3-6] or persistent ground state 
currents in superconducting rings [7, 8]. 

Several interesting effects in condensed matter experiments 
typically require a system to be exposed to a magnetic field. 
Therefore, a quantum simulator of these effects needs to in- 
clude or simulate such a magnetic field. Ultracold atoms offer 
several promising ways to achieve this goal. One idea is to 
rotate a gas of atoms close to the frequency of the harmonic 
confinement that traps the atoms. This may result in states 
also found in the quantum Hall effect [9-11]. However, the 
high rotation frequency ft of the atoms introduces a large cen- 
trifugal term, which tends to destroy the atomic cloud. On the 
other hand, to reach the strongly correlated state the rotation 
frequency of the cloud and the frequency of the confining har- 
monic trap ujc have to be balanced with high accuracy such 
that 1 - Z/2Nas < il/uJc < 1 - Shs/ZN, where Z is the 
transverse distance over which the condensate is uniform, N 
is the number of atoms, and is on the order of the s-wave 
scattering length [12]. The necessity of precise counterbalanc- 
ing is alleviated, in principle, by using traps with potentials 
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Stronger than harmonic at the expense of harder experimen- 
tal implementation [13]. A different approach is to use lasers 
with an orbital angular momentum [14, 15] or to use Raman- 
assisted hopping in a layered optical lattice geometry [16] in 
order to mimic the effect of a magnetic field. Both proposals 
might also be difficult to realize experimentally. 

In this paper, we analyze and extend methods introduced 
by some of the present authors [17] for simulating a mag- 
netic field in atoms trapped in a ring-shaped, stationary op- 
tical lattice. These lattice atoms, in the following also 
called impurities, are immersed into a ring-shaped quasi- 
onedimensional Bose-Einstein condensate (BEC) exhibiting 
a nontrivial phase. For creating this phase we consider two 
possible setups. Firstly, a BEC trapped in a rotating ring [18- 
20]; secondly, a BEC carrying a quantized angular momentum 
in a stationary ring [21], i.e., a BEC in an excited state. In the 
former case, we can allow a confining potential to overcom- 
pensate the centrifugal term in the rotating BEC as we are 
mainly interested in the dynamics of the impurities. This is in 
contrast to methods for observing quantum Hall states directly 
in a rotating gas of ultracold atoms where the centrifugal term 
needs to be counterbalanced precisely [12]. The BEC is de- 
scribed in terms of Bogoliubov excitations, which couple to 
the impurity atoms. It has been shown in earher publications 
that the presence of the phonon-impurity coupling leads to ob- 
servable effects in the transport properties of the lattice atoms 
such as the transition from coherent to incoherent transport 
with increasing temperature or clustering [22-24]. In addi- 
tion, by allowing the BEC to exhibit a nontrivial phase the 
coupling leads to an artificial magnetic field in the dynamics 
of the impurities which manifests itself as a nontrivial phase 
term in the hopping of the lattice atoms. By describing the 
impurities in terms of polarons we derive an analytical ex- 
pression for the induced phase. 

The induced phase changes the dynamics of the lattice 
atoms considerably. We show that it leads to current carry- 
ing ground states in the impurities, which can be detected ex- 



2 



perimentally. Another experimentally accessible method for 
detecting the phase of ultracold atoms is to evaluate time of 
flight (ToF) images [25]. We discuss this method as a means 
of revealing the rotational state of the underlying BEC with 
the impurity as a probe. We find that nonzero temperature and 
increasing interspecies coupling obscure this transition. We 
also show that the tendency of the lattice atoms to assemble in 
clusters depends on the value of the induced phase. Further- 
more, we observe the onset of quantum chaos if the symmetry 
of our two-species system is broken by a disorder potential. 
Quantum chaos is characterized by the distribution of energy 
levels. If this distribution is of the Gaussian orthogonal en- 
semble (GOE) type, then the system is said to behave quan- 
tum chaotically [26, 27]. Although we restrict ourselves to a 
quasi-onedimensional system throughout this paper, our oned- 
imensional (ID) model will allow insights into experimentally 
more realistic twodimensional (2D) setups. 

This paper is organized as follows. In Sec. II we derive 
an effective Hamiltonian for the impurity atoms immersed 
into the BEC. The impurities in this model are dressed with 
a cloud of BEC excitations, which induces additional in- 
teractions and a phase term on the impurity hopping. The 
dressed impurities can be interpreted as polarons [22, 28]. 
We derive analytical formulas for the induced phase in dif- 
ferent experimental setups. In Sec. Ill we present two meth- 
ods for detecting the additional phase twist in the impurities. 
These methods are based on ground state current measure- 
ments and ToF imaging. In Sec. IV we study the effect of the 
phonon-mediated interaction in the impurities on their clus- 
tering properties. Finally, we extend previous studies [29] of 
the chaotic behavior of a one-species Bose-Hubbard model 
with a phase twist to two species in Sec. V. We conclude in 
Sec. VI. The application of the Bogoliubov approximation to 
the Bose-Hubbard model with a phase twist is discussed in an 
Appendix. 



II. THE MODEL 

We consider two species of ultracold atoms, species a and b, 
both confined in a quasi- ID ring geometry. Atoms of species 
b form a BEC confined in a ring trap. Atoms of species a, 
considered as impurities in the BEC, are trapped in a ring- 
shaped optical lattice [30]. This optical lattice is submerged 
into the BEC which is not affected by the presence of the 
optical lattice potential [31]. For a detailed derivation of the 
model in higher dimensions we refer to Refs. [17, 22, 23]. The 
term quasi- ID is taken to mean that the atoms are confined in 
a tight transverse trap such that they are in the ground state 
of this trap. A transverse harmonic trap with frequency uj± 
introduces the characteristic length scale £± = \JhJmy^, 
where m\y is the mass of atoms of species b. For the gas to 
be in the transverse ground state of this trap we require that 
huj\_ ^ gn, where g = 2lT?as/m\jl\ is the ID interaction 
strength with the three-dimensional (3D) s-wave scatter- 
ing length and n = N/L the linear density of the gas in 
a ring of length L with N atoms. Here we have assumed 
that the scattering of the atoms is still essentially 3D, i.e.. 



as. In lower dimensional gases phase and density fluc- 
tuations play an important role. Our calculations are valid 
if the phase of the BEC does not fluctuate over the circum- 
ference of the ring L. In a ID BEC coherence is preserved 
over a length scale — nfi^ /iribkBT, where /cs is Boltz- 
mann's constant and T is the BEC temperature [32]. As we 
will employ the Bogoliubov approximation, density fluctua- 
tions must be small. This is ensured in the weakly interacting 
limit n^/i ^ 1, with — f^l ^Jrubgn the healing length, if 

T <t^Td = ^K^n^g/rrni/kB [33]. In this temperature regime 
the atoms are in a BEC because Td lies below the condensa- 
tion temperature [34, 35]. 

We write the Hamiltonian of the full system as H = Hb + 
Ha+Hi, where the three parts represent the BEC, the impuri- 
ties, and the interaction between the two species, respectively. 
The BEC Hamiltonian Hb and interaction Hj are given by 

Hb = Jdx^^x) [Hq + ^^\x)^{x)] ^(x), 
Hi = K dxx'{x)x{x)4'''{x)^{x). 



The field operator 4>{x) annihilates a BEC atom at position x, 
whereas x{^) annihilates an impurity at x. The Hamiltonian 
Hq describes the noninteracting part of the BEC and will be 
defined when we discuss specific examples, g is the ID inter- 
action strength of the BEC atoms and k is the ID interaction 
strength between a BEC atom and an impurity. We assume 
that the temperature is sufficiently low such that the impurities 
in the lattice only occupy the lowest band. Then we can ex- 
pand their field operator in terms of Wannier functions rjj {x) 
as [36] 
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The operator aj annihilates an impurity in lattice site j. We 
assume the lattice to be sufficiently deep such that the tight- 
binding approximation for the Wannier functions holds, i.e., 
j dx\rjj{x)\'^\T]j'{x)\'^ ~ for different lattice sites j ^ j'. 
With these assumptions the impurities are well described by a 
Bose-Hubbard Hamiltonian [37] 



(3,3') 



3 



where Jq is the hopping matrix element between two neigh- 



boring sites, Ua the onsite interaction, hj 
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and fla 



is the chemical potential. We note that it is straightforward 
to extend our derivation to fermionic impurities by using a 
Fermi-Hubbard Hamiltonian. 

We write the BEC field operator as (t){x) = 4>o{x) + C{x), 
where (j)o{x) is the solution of the Gross-Pitaevskii equation 
(GPE) for K = 



[Ho + g\M^)\^]M^) = 



(1) 



and ({x) describes a small quantized deviation from the 
mean field solution (f>o{x). We will be studying rotating 
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BECs, which means that 4>o{x) is, in general, a complex 
function. Furthermore, we assume that the impurity-boson 
coupling K fulfills \K\/gno{x)£,h{x) ^ 1, where no(x) — 
|(/)o(x)p is the mean field density and we allow a space- 
dependent density in the healing length. The smallness con- 
dition for the interspecies coupling k ensures that the defor- 
mation of the BEC due to the presence of the impurities is 
small [23]. We plug the expansion 0o(a:) + ({x) into the 
Hamiltonian Hb + Hj and keep terms up to second order 
in K. The linear terms in ({x) result in the interspecies in- 
teraction K J dxX' {x)x{x)[(j)Q{x)(^ {x) + (I)q{x)(^{x)]. The 
quadratic terms in (^{x) are diagonalized by expanding the 
fluctuations in terms of Bogoliubov excitations as ({x) — 
[uq{x)bq — v*{x)b']j]. Here bq annihilates an excitation 
(phonon) in mode q and the prime indicates that the sum does 
not include the BEC mode. The operators bq and 6j fulfill 
the usual bosonic commutation relations [^^,6^/] — Sqq> and 

[bq, bq'] = 0. The mode functions Uq{x) and Vq{x) satisfy the 
Bogoliubov-de Gennes (BdG) equations [17] 



allows us to write down the effective Hamiltonian 

H^ff = -Ja E e'""-^'a]a,, + ^^^^^^^ " 1) 
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The chemical potential is jij = jj.a — ktio 



(4) 



Uj = 



Ja exp - 



Ua — "^Ej, the reduced hopping J a 



'[2NqiT) + l]/2 



where 



Nq{T) = [exp(hujq/kBT) - and the interaction 

Vjjv = J2' hiOq{Mj^qM*,q + M*qMf^q). The Boltzmann 

distribution Nq{T) enters because we have averaged the 
BEC degrees of freedom over a thermal phonon distribution 
at temperature T. The phase factor in the hopping term of 
Eq. (4) is given by 
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[Ho + 2gno(x)]uq{x) - g[(l)o{x)]'^Vq{x) = hwqUq{x), (2a) 

[hI + 2gna{x)]vq{x) - g[(j)o{x)]'^Uq{x) = -hWqVq{x). 

(2b) 



The eigenenergies of these equations are the quasiparticle 
energies huq. The total Hamiltonian is then given by the 
Hubbard-Holstein Hamiltonian [38] 



H = Ha + J2Yl' f^ii^l3A + Mlqbl)7 

3 1 

+ Y^jnj + ^ hujqblbq. 



(3) 



Here, Mj^q — {n/hiUq) J dx[(j)Q{x)uq{x) — 
4'o{x)vq{x)]\'i]j{x)\'^ are matrix elements of the phonon- 
impurity coupling, and Ej — k J dxno{x)\rij{x)\'^ is a mean 
field shift. By using the Lang-Firsov transformation [28] 
this Hamiltonian can be shown to describe the dynamics 
of polarons in an optical lattice [17, 22]. The polarons 
in this model are given by the impurity atoms surrounded 
by a coherent cloud of Bogoliubov phonons. In the fol- 
lowing we require that fcsT <C Ej and Ja <C Ej, where 

El 
hujq\Mj^q\'^ is the polaronic level shift. For 

the cases considered in this work it is independent of the lat- 
tice site j. We assume the phonons to be thermally occupied 
at temperature T. Furthermore, the characteristic hopping 
speed has to fulfill dJa/h ^ c, where c ^ \/gno/mb is the 
speed of sound and d the lattice spacing. These conditions 
ensure that the polaron dynamics is coherent, retardation 
effects are suppressed, and that the hopping term can be 
treated as a perturbation [22]. As shown in Ref. [17] this 



Since we allow the condensate wave function 4'q{x) to be 
complex, in general, ajji is nonzero. This means that the 
impurity atoms pick up a phase when they hop across lattice 
sites. The derivation above demonstrates that the origin of this 
induced phase is a coupling of quasiparticles in the BEC to 
the impurities. The occurrence of such a phase in the Bose- 
Hubbard model leads to a change in the dynamics of the impu- 
rities in the lattice. In the following we will consider concrete 
systems and derive the analytical expressions for the corre- 
sponding induced phases. 



A. The BEC in a rotating ring 

For the first model system we restate here the results of [17] 
for a BEC in a rotating ring and also introduce quantities nec- 
essary for the discussion in the remainder of the paper The 
quasi- ID BEC is confined to a ring of length L rotating at 
angular frequency fl. The impurities are trapped in a ring- 
shaped optical lattice, which is submerged into the BEC. The 
noninteracting part of the Hamiltonian of the BEC is 



2mb dx'^ dx 



rot 



(6) 



where v = Rft is the rotational speed, R = _L/27r is the ra- 
dius of the ring, and the chemical potential of the BEC. 
The function (I)q{x) = y/no exp{ikx) solves the correspond- 
ing GPE, Eq. (1), with the quantized momentum k = 2-Km/L 
and m an integer such that k = mi^v/h — Ak, where Afc e 
[— tt/L, tt/L). The density of the BEC is denoted with uq. 
This definition ensures that the BEC is in the ground state 
with its momentum closest to the momentum of the rota- 
tion m}jv/h. The chemical potential of the BEC is given 
by Mb" = fi.^fc^/2mb + griQ — hvk. The fluctuations around 
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the ground state wave function lead to Bogoliubov excitations 
with energy 



fv^qAk 

nib 



(7) 



where = ^e°(e° + 2gno) and = h'^q^/2mi,. This Bo- 
goliubov phonon energy contains the mismatch Afc between 
the ground state momentum and the angular momentum given 
by the rotation speed. We rewrite the definition of the mis- 
match as Aw := hAk/mi, = v — hk/rrn, —: v — Vb, where Vb 
is the velocity of the EEC. Positive excitation energies require 
that lAul < miug^o /hq ^ c. This asymptotic behavior of 
the velocity follows if we only consider low lying excitations 
(since L ^ The Bogoliubov modes satisfy the periodic 
boundary conditions, which results in the quantization of their 
quasimomenta as g = 27Tm/L with m an integer. The solu- 
tions of the BdG equations, Eqs. (2), are the mode functions 



Uq{x) 



(8a) 
(8b) 



where = 1 - 
cients fulfill Ug ± Vq 



= [(eO + gno)lE^ + l]/2. Their coeffi- 
(i?^/e")='=^/^. In a tighdy confining 
lattice the Wannier functions of the impurities are well de- 
scribed by Gaussians rij{x) = exp[— (x — XjY /2a'^]/ yj -jTia 
with width a centered at a lattice site j. The position of the 
j-th lattice site is parametrized by its angle on the ring, 
i.e., Xj = ^jNsd/2'K, where Ng is the number of lattice sites 
and d is the lattice spacing. The coupling matrix elements are 
then given by 




iqxj-q 



V4 



The Gaussian width a of the localized impurities introduces 
a cutoff in the phonon momenta contributing to the coupling. 
By plugging this result into Eq. (5) we finally arrive at the 
induced phase of the impurities 
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27rL ^ EB [hiofY 



sin((7d), (9) 



where a™' = aj+ij since we only consider nearest-neighbor 
hopping. As a result of the momentum cutoff in the coupling 
terms, the momenta contributing to the sum in the phase are 
distributed with a Gaussian envelope. The reduced hopping 
term is given by 



In Fig. 1(a) we plot the induced phase for typical ex- 
perimental parameters for different BEC coupling strengths 
g [39, 40]. The values of the induced phase a™' ~ 0.03 are 
sufficiently large to cause observable effects in the dynamics 
of the system. In our coherent approximation, the effect of 
nonzero temperature is only to reduce the effective hopping 
rate. For simplicity we have set the temperature to zero in the 
figure and the magnitude of the hopping J"' does not deviate 
significantly from the magnitude of the bare hopping J^. In- 
coherent hopping will only affect the phase if the temperature 
lies above the polaron energy [22]. Density fluctuations will 
play a significant role for temperatures T > Td- If the length 
of the ring exceeds the coherence length Ic, the assumption of 
a well defined phase over the whole BEC is invalidated. The 
temperature scale for thermal phonon excitation is of similar 
magnitude. Furthermore, temperatures at fc^T > mrh^ /nibR 
will excite additional rotations of the BEC, which results in 
a multi-peaked distribution of rotation velocities and a grad- 
ual loss of superfluidity [41]. Typically, all these effects start 
playing a significant role at temperatures on the order of a few 
tens of nanokelvin. Hence, the effects discussed here can be 
observed for sufficiently cold samples. 

From Fig. 1(a) we see that at first the induced phase de- 
pends linearly on the rotation frequency of the BEC. At the 
critical frequency ilcrit = h/2mbR^ we observe a sudden 
jump to negative phases. This jump is caused by a jump in 
the momentum of the BEC at the critical rotation frequency. 
Because of the quantization of the quasimomentum in units of 
2tt/L, at this frequency the ground state of the BEC changes 
to a nonzero quasimomentum. These jumps repeat at odd mul- 
tiples of the critical rotation frequency. 



B. The BEC in an optical lattice 

In order to be able to compare our analytical results with nu- 
merical solutions we now introduce a system which is numer- 
ically solvable beyond the Bogoliubov approximation. This 
allows us to study the validity of our model of an induced 
phase beyond the Bogoliubov approximation. We assume that 
the BEC consisting of N^. atoms is trapped in a ring-shaped 
optical lattice with Ng lattice sites and spacing d, which is 
independent of the lattice for the impurities. We will derive 
an analytical expression for the induced phase in the Bogoli- 
ubov approximation, which will have a similar form as the 
expression in the continuous case of Sec. II A, Eq. (9). We 
will use these results for the discussion of the numerical re- 
sults in Sec. III. One can show that in the limit of A^^ — > oo 
with Ngd = L = const the two expressions coincide. 

The Hamiltonian of the condensate is given by a general- 
ized Bose-Hubbard Hamiltonian [42] 
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(10) 



Uc 



(11) 



5 



(a) 



n (rad/s) 
100 200 



0.02 



0.00 



-0.02 
(b) 0.02 



-0.0 




0.02 



0.04 



0.06 



FIG. 1: (Color online) Induced phase in the lattice with 30 lattice 
sites for (a) a ring rotating with angular speed Q and (b) the BEC 
in a lattice with phase term Qc. We assume a BEC of *^Rb with a 
linear density no — 5 x 10^ m~^. The impurities are *^Rb atoms 
in a different hyperfine level. The lattice spacing is d = 400 nm, 
which means that the ring has a circumference of L = 12 fim. In 
(a) the interspecies coupling strength is K,/2dER — 0.045, where 
Er — h^ix"^ /2mad? is the recoil energy of the impurity lattice with 
impurity mass rUa, and the curves indicate different BEC couplings 
g/2dER = 0.01,0.015,0.02 (solid, dashed, dashed-dotted lines, 
respectively). In (b) we assume that both lattices have the same 
depth, interspecies coupling is Ui/ Jc = 1-5, and BEC couplings 
Uc/ Jc = 0.5, 1, 1.5 (solid, dashed, and dashed-dotted lines, respec- 
tively). 



where Jc and Uc are the hopping and onsite interaction 
terms of the BEC, respectively, and Cj is the bosonic anni- 
hilation operator for site j. For a lattice rotating with an- 
gular frequency VI the phase factor ac is given by Uc — 
mb^Nsd? I ^Tr'^h [42]. The impurities are trapped in a second 
lattice as before, for which we assume that it has the same 
spacing d and the same number of lattice sites as the BEC 
lattice. In this picture the interaction between impurity and 
BEC atoms is given by 



(12) 



where Uj is the coupling strength between the two species. In 
the remainder of this section we will derive analytically the 
induced phase of this system resulting from our approximate 
polaron model. Later we will use Eqs. (11) and (12) for nu- 
merical computations of the full two-species model. 

Similar to the implementation discussed in the preceding 
section we now expand the annihilation operators Cj around an 



order parameter as Cj 



+ ^j. Then we rewrite this ex- 



pansion in the momentum representation with quasimomen- 
tum 2'nq/Ng, where q is an integer, and apply a Bogoliubov 
transformation with quasiparticle operators hq and For the 
fluctuations to be small we require that Uj /UcfiQ <^ 1, where 
fio = Nc/Ns is the number of BEC atoms per lattice site. 
The main result is that, as before, we aiTive at an effective 
Hamiltonian of the Hubbard-Holstein form with polarons as 
the quasiparticles hopping across lattice sites. A derivation of 
the results is given in the Appendix. In the coherent dynamics 
limit the Bogoliubov excitations mediate additional interac- 
tions and induce a phase in the impurities. The excitations in 
this system have an energy 



huji 



Ei 



A„ 



(13) 



where E„ 



and 



-2 Jc sin{2Trq/Ns) sin(A6'/7Vs), 
4JcSin^(7rg/iVs)cos(A6l/A^s). The angle Ae* is deter- 
mined via the relation 27rac = 2ttv/Ns + A9/Ns with h' 
an integer and Ad E [— 7r,7r). The ground state momentum 
is determined by the winding number v = [acNs] , where 
the symbol [•] denotes rounding to the nearest integer. One 
can show that the phonon energies in Eq. (13), reduce to the 
excitation energies of the rotating ring, Eq. (7), in the limit 
Ns —i' oo with Ngd — L — const. The coupling matrix 
elements between impurities and phonons are given by 



M. 




g_^27viq]/N, 



(14) 



From Eq. (5) we calculate the induced phase of the impurity 
atoms as 



tt2- ^=-1 

9=1 



,L\2 



sm{2Trq/N,). 



(15) 



The form of this induced phase is very similar to the one for 
a BEC in a rotating ring, Eq. (9). Here however, the cutoff in 
the phonon quasimomentum is introduced by the quantization 
of the momentum through the finite lattice. This results in a 
finite sum instead of a Gaussian depending on the impurity 
locaUzation width a. The reduced hopping term is given by 



Ja = Ja exp 



q=l 



1 



E^ {ha 



.L\2 



X [l-cos(2^g/iV,)][2iV,(r) + l] 



Numerical tests show that the rescaled induced phase 
a^27TNs/U]no becomes maximal for A9 — >■ tt and C/c^o — >■ 
0. However, Bogoliubov theory loses validity at A9 — tt be- 
cause the ground state becomes degenerate. At small Ud Jc 
we find Tlq :^ n = Nc/Ng so in order to maximize the con- 
stant factor in the induce phase, we have to choose a large 
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operator is then given by 
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FIG. 2: (Color online) Induced phase with a BEC in an excited state. 
The parameters are the same as in Fig. 1(a) and we abbreviate g = 
g/2dER. The critical momenta for the different interaction strengths 
are given by hn,L/2-K ~ 2.1, 2.6, 3.0 (for g = 0.01,0.015,0.02, 
respectively). The lines are only to guide the eye. 



filling n and large interspecies coupling [//. However, the in- 
terspecies coupling has to remain sufficiently small to fulfill 
the condition Uj/UcTi <C 1 that was necessary for the analyt- 
ical derivation. 

In Fig. 1(b) we plot the induced phase according to Eq. (15) 
for similar parameters as in Fig. 1(a) but now with a BEC in 
a rotating optical lattice. The induced phase now critically 
depends on the BEC phase ac, which takes the role of the 
frequency in the preceding section. We observe jumps in the 
induced phase at critical BEC phases, which have the same 
origin as the critical frequencies in the rotating ring. The crit- 
ical phases occur at acnt = (2j -f l)/2iVs, where j is an in- 
teger. We can expect similar orders of magnitude for the in- 
duced phase for both setups with or without an optical lattice 
for the BEC atoms. Note that both the upper and lower x-axes 
in the figure are valid for both plots because we assume two 
species of the same mass and the phase in the Bose-Hubbard 
Hamiltonian Hbh is a linear function of the frequency. 



C. The rotating BEC in a static ring 

Again we consider a quasi-lD BEC in a ring of radius 
R = L/2tt and the impurities trapped in a ring-shaped op- 
tical lattice immersed into the BEC. In contrast to Sec. II A, 
here we assume that the BEC rotates but the ring does not. 
Hence, the BEC is not in the ground state of the Hamilto- 
nian. This excited state can be created, for example, with a 
STIRAP process [15, 21] or by initializing the BEC in a ro- 
tating ring whose rotation is then turned off for the duration 
of the experiment. The interaction-free Hamiltonian has the 
same form as the one in the rotating ring, Eq. (6), but with 
?; = 0. As before we expand the fluctuations around the order 
parameter in terms of Bogoliubov excitations. The BEC field 



Akx 



(16) 



Periodicity requires the quantization of the phase k = 2-Km/L 
with m an integer. Plugging this ansatz into the GPE results 
in the chemical potential /i^""^ ~ h?k'^ /2mi, + guQ. The so- 
lutions of the BdG equations diagonalize the quadratic pait of 
the Hamiltonian. In [43] they are given as 



Uq{x) 



Vq 



with the same coefficients Uq, Vq as in the case of a rotating 
ring [see definition below Eqs. (8)]. The energies of these 
solutions are given by 

huj'g'"' = E^ + hqw, 

where w — hk/rrn, is the speed associated with momentum k. 
For fiij^"'^ to be positive we have assumed that the velocity of 
the BEC is less than the speed of sound, i.e., w < c, which 
is equivalent to a critical momentum fcd-it = Again we 

assume that the impuiity wave function is given by a Gaussian 
centered at the position Xj of lattice site j. This allows us to 
evaluate the induced phase as 



2ttL ^ 

q 



1 



sinM). 



We see that the form of this phase is the same as in Eq. (9) but 
now the phonon energy depends on the full momentum of the 
BEC k, not only the mismatch to the ground state Afc. Simi- 
larly, the form of the reduced hopping is the same as Eq. (10) 
but with /lo;"' replaced by /iw^"'^. The higher energies allow 
for larger induced phases but entail the experimental difficulty 
of maintaining the BEC in an excited state [8, 44]. A deriva- 
tion of the life time of this excited state is beyond the scope of 
the present work. 

Figure 2 shows the induced phase in this setup for experi- 
mentally accessible parameters below the critical momentum. 
Induced phases up to a™'^ ~ —0.25 are readily achievable 
in this setup and even higher phase are possible by slightly 
changing the parameters. The reason for the negative sign 
of the phase is that the dispersion favors phonons with nega- 
tive quasimomentum. A negative quasimomentum, i.e., total 
phonon quasimomentum k — \q\, pushes the system closer to 
the ground state at zero total momentum. 



III. PROBING THE PHASE TWIST 

In [17] it was shown that the presence of a phase twist in 
the impurity Hamiltonian can be observed as a drift of the 
impurities. If we prepare a Gaussian-shaped distribution of 
impurity atoms centered at lattice site j, then it will expand 
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and the mean position of this packet will drift either to the left 
or right depending on the induced phase. This drift can be 
verified experimentally. It is also present for impurity speeds 
below the speed of sound, which, at first glance, seems to con- 
tradict the Landau criterion of superfluid flow [45, 46]. From 
this criterion one would expect that a BEC rotating below the 
Landau critical speed does not impart a momentum on the im- 
purity. Recent work by Sykes et al. [47] and earlier works by 
Roberts et al. [48, 49] suggest Doppler-shifted scattering pro- 
cesses of quantum fluctuations in the BEC with the impurity 
as a source for this drift. When an incoming wave of a BEC 
quantum fluctuation is reflected off an impurity it experiences 
a Doppler shift depending on the propagation direction, which 
leads to an overall drag force on the impurity. In contrast to 
their derivation, our results are based on the nonperturbative 
Lang-Firsov transformation and we assume periodic boundary 
conditions. 

In the remainder of this section, we present detection meth- 
ods for the induced phase twist, namely, via mass current 
measurements and ToF expansion. In order to compare the 
effective one-species model, Eq. (4), with numerical results 
we use a two-species Bose-Hubbard Hamiltonian H2BH = 
Ha + Hbh, which can be solved by numerical diagonaliza- 
tion for small systems. 



A. Impurity ground state momentum 

For a single impurity our effective one-species model pre- 
dicts a vanishing ground state momentum for la^l < l/2Ng. 
In the full two-species model interactions between the two 
atom species cause a broadening of the reduced impurity 
ground state in momentum space even for a single impu- 
rity. We now investigate how this broadening can be used 
to measure the induced phase aa of our effective model. For 
vanishing induced phase the broadening in momentum space 
is symmetric around the zero momentum state. However, 
for < \aa\ < l/2Ng the broadening becomes asymmet- 
ric around the zero momentum state. This leads to a small 
nonvanishing mean momentum of the impurity ground state. 
We note that the origin of this effect is different from the 
expected jump of the momentum ground state at the criti- 
cal phases aci-it- At a critical phase the macroscopic occu- 
pation of a momentum state changes, whereas the asymmet- 
ric broadening is a perturbative effect and does not cause 
such a macroscopic shift. To see this effect analytically we 
calculate the ground state of the Hubbard-Holstein Hamilto- 
nian, Eq. (3), in first-order Rayleigh-Schrodinger perturba- 
tion theory for a single impurity [50]. For simplicity we let 
Ej = and /ia = so that we can rewrite Eq. (3) in mo- 
mentum space as iJnH = Z]fe -^fe^l^fc + X/ ^^THK + 



k+q^k 



where dk an- 
nihilates an impurity with quasimomentum k, E'j, = 
— 2 Jq cos(fcc?) is the unperturbed impurity dispersion with 
quasimomentum k, and il/"* is the Fourier-transformed cou- 
pling matrix element Mj,q. The coupling term mixes phonon 



and impurity excitations. We write the unperturbed states as 
\k,Nq), which indicates an impurity with quasimomentum k 
and Nq Bogoliubov excitations with quasimomentum q. With 
this notation the perturbed (dressed) ground state is 



|/=o,0) 



(1) _ 



hjJqMq 



hui, 



\ko - q, Iq) , 



where ko is the quasimomentum of the impurity ground state. 
The quasimomentum distribution and mean quasimomen- 
tum k of this state are given by 



/~t- \ X f f>uJko-k\Mk„^k\ 



k = y^nfcA;. 



, (17) 



(18) 



To be able to compare the numerically obtained mean quasi- 
momentum with this formula we insert the expression for Mk 
for the BEC in an optical lattice, Eq. (14), into Eq. (17). 
Moreover, in a lattice the quasimomentum k is expressed 
as 2TTk/Ns, where k is now an integer. With the resulting 
quasimomentum distribution nt the mean quasimomentum, 
Eq. (18), reads 



Ufno 



fc=-LW,/2]+l 



2TTk/Ns 



Ek (E^ 



(19) 



where we have used the fact that e^, Ej^, and E^ are even 
functions of k. This expression should be compared to the 
expression for the induced phase, Eq. (15). The sine in 
Eq. (15) can be linearized for 2TTq/Ns close the roots of the 
sine. Furthermore, the summation indices of both expressions 
can be shifted so that the two expressions formally coincide 
apart from the different energy denominators and a constant. 
To analyze the energy denominator in Eq. (19) we assume 
Uc = for clarity. Then it is given by Eq — E^ — huj^ = 

-2{Ja + Jc)[l - C0s(2^A:/iV,)] - Afe ~ -/i^^ll + Ja/Jc) 

for phases close to j/Ng (j integer), that is small Aj,. We 
might thus expect that 



27r 



(1 + Ja/JcY 



(20) 



In order to check the validity of these approximations 
we have computed ground states of the two-species Bose- 
Hubbard Hamiltonian with a single impurity. From the re- 
duced density matrix of the impurity we compute the mean 
quasimomentum and compare it with Eq. (19). Figure 3(a) 
shows the error of the formula (19) relative to the numerical 
results for different sets of interactions. As expected, the ac- 
curacy of the analytical result is better for small interactions 
(e.g., approximately ±5% for Uc — O.Uc)- In Fig. 3(b)- 
(c) we compare the mean quasimomentum of the impurity 
ground state with the analytical result for the induced phase, 
Eq. (15). Clearly, our approximations discussed in the pre- 
ceding paragraph predict the value of the induced phase for 
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FIG. 3: (Color online) (a) Relative error /Cana / ^niim 

1 of 

Eq. (19) compared to the numerically exact ground state quasi- 
momentum for different interactions at ~ 12, 6 BEC 

atoms, and one impurity with Ja ~ Jc- Line styles indi- 
cate different BEC interactions Uc (see legend) and increasing 
lightness (from top to bottom) indicates decreasing inter- 
species coupling Ui/Jc e {0.2,0.15,0.14, ...,0.05,0.01} 
(dashed-dotted lines), {0.9, 0.8, . . . , 0.1, 0.05, 0.01} (dashed), 
{1,0.7,0.6, ... ,0.1,0.05,0.01} (solid), (b)-(c) Comparison of 
analytical induced phase (solid), numerically exact mean quasimo- 
mentum (dashed) and the rescaled analytical mean quasimomentum 
-2fc/7r (dashed-dotted) for (b) Uc = 0.5Jc, Ui = 0.4 Jc and (c) 
Uc = Jc, Ui = 0.2Jc. The rescaling factor —2/ii results from 
Eq. (20). The dotted vertical line indicates Uc = l/Ns- 



close to Oc — j/Ng (j integer), that is away from the crit- 
ical phases a„it = (2j 4- l)/2iVs, where Bogoliubov theory 
fails. Also note that in Fig. 3(b) the constraints on the validity 
of our analytical derivation of the induced phase are not met 
(Ui <4c Ucfio)- Nevertheless, the rescaled mean quasimomen- 
tum approximates the induced phase well close to the roots 
of a^, which means that qualitatively our effective polaron 
model still describes the underlying physics. 

Our computational resources restrict the size of the sys- 
tems studied here to moderately small number of lattice sites 
and filling factors. Therefore, we could not directly observe 
the jump of the impurity to a macroscopic occupation of a 
nonzero momentum ground state for > Qciit- However, 
our results clearly indicate an imbalanced momentum distri- 
bution owing to the influence of an induced phase. For larger 
number of lattice sites the critical phase a„ii decreases. Fur- 
thermore, the quantization of the impurity ground state mo- 



mentum in units of 2tt/Ns becomes denser for larger lattices 
so that more allowed momentum states become available for 
the impurity. Therefore, we expect to see a macroscopic pop- 
ulation of a nonzero momentum state in larger systems caused 
by a large induced phase as predicted by our effective model. 
Such a macroscopic population will be observable as an im- 
purity current in an experiment. The current measured by the 
operator 



Ca 



iJama 



E 



is closely related to the mean momentum of the impurity. It is 
straightforward to show that the cuiTent of a single impurity 
with momentum k is (Co) = 2 Jq ma sin(27rfc)//i. For small 
k as in Fig. 3 we may approximate (Co) oc sin(27rfc) ~ 27rfc. 
In the preceding paragraph we have seen that close to the zero 
crossing the induced phase and mean momentum are propor- 
tional, which shows that the induced phase and the mass cur- 
rent measured by Ca are proportional around = l/iVg. 



B. Time of flight images 

Another experimentally readily accessible method to probe 
for a phase term is the ToF expansion of the atoms. After their 
evolution the atoms are abruptly released from the trap and 
expand for a time t before they are imaged. If interactions 
can be neglected during the time of flight, the imaged density 
distribution coiTesponds to the momentum distribution of the 
atoms in the trap. This distribution at a position x is given by 



Ptof(x) = \w(K.] 



jl J2 



where we neglect a constant factor {ma/ht)^ [51]. In the bal- 
listic approximation the momentum is K = niax/ht with 
the mass of an impurity atom, and 7«(K) the Fourier trans- 
form of the Wannier function of the impurity trapped in an 
optical lattice. The one-particle density matrix is given by 



/9i(Rjj,Rj2) = {a^. Or. ), where a^, creates an impurity 
atom at position R^ . In a ring of radius R the positions Rj 
are fixed at R^ = R[sm(2Trj / Ng) , cos{2Trj / Ng)]. For tightly 
localized atoms we can set w(K) = l/^/Wl and write 



N. 

X exp{iis:^i?[sin(27rji/iV,) - sin(27rj2/A^.)] 
+ iKyR[cos{2nji/N,) ~ cos(27rj2/iV,)]}. 

To illustrate the effect of an induced phase term in the ef- 
fective impurity Hamiltonian, Eq. (4), we first consider the 
ToF expansion of a single atom, which is not immersed in 
a BEC. While this atom alone without a suiTounding BEC 
would not exhibit a phase term in practice, it is possible to 
include this phase term in the numerical calculation to illumi- 
nate its effect. We will study the full two-species system in the 
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FIG. 4: (Color online) Time of flight distribution pToF of a single 
impurity released from a ring lattice with A^^ = 12 at T = in 
the single-species model, Eq. (4). The central peak in (a) indicates 
a zero quasimomentum ground state for Ua ~ 0.04. In (b) the van- 
ishing density in the center indicates that the ground state exhibits 
nonzero momentum {Ua = 0.05). The jump between these two mo- 
mentum states occurs at aa = l/2Ns — 0.042. In (c) we show the 
corresponding density profiles indicated with horizontal lines in (a) 
and (b). The ToF expansion time is t = 50 ms for ^'^Rb atoms in a 
ring of radius R — 12/im. 



next paragraph. We have seen in the preceding sections that 
the critical phases acrit correspond to a macroscopic jump in 
the momentum of the ground state. As a consequence, if the 
induced phase crosses aan, the ground state of the impurity 
will jump to a higher momentum state. Since the ToF dis- 
tribution represents the momentum distribution of the trapped 
atoms, we expect a central peak for nonrotating impurities and 
a vanishing density at the center for impurities with a nonzero 
momentum [25, 52]. ToF images can reveal this jump in the 
momentum of the impurity with good accuracy. Correspond- 
ing density plots are shown in Fig. 4(a) for a phase below 
and in Fig. 4(b) for a phase above the lowest critical phase 
ttciit = l/2iVs. In Fig. 4(c) we plot the profiles of the density 
distribution after a realistic time of flight. The different pro- 
files of the two ToF images should be clearly distinguishable 
with typical camera resolutions [53]. 

We now return to the numerical study of the full two- 
species Bose-Hubbard model with a rotating BEC and a sta- 
tionary impurity lattice. For the moderately small system 
sizes studied here we cannot expect to see a clear transition 
to a macroscopic occupation of a nonzero impurity momen- 
tum state because the induced phase is not sufficiently large 
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FIG. 5: (Color online) Influence of temperature on the ToF images 
in the two-species Bose-Hubbard model. A single impurity is sub- 
merged into a BEC with 6 atoms in Ns — 12 lattice sites, hop- 
ping Ja = Jc, interspecies coupling Ui — 0.5Jc, BEC interaction 
Uc = OAJc, and BEC phase = 0.042. In (a) kgT = and in 
(b) fcsT = 2Jc. 



to cause this jump. Instead, in Fig. 5 we show the influence 
of nonzero temperature on the central peak in the ToF image 
of the impurity. The BEC phase is chosen to be close to the 
critical phase 1/2 A^^ so that the induced phase is expected to 
be large. The density distribution in Fig. 5(a) is very simi- 
lar to the single-species calculation in Fig. 4(a), which means 
that small interspecies interactions do not influence the ToF 
image significantly. However, the induced phase is not suffi- 
ciently large to cause a macroscopic occupation of a nonzero 
momentum state as revealed by a vanishing central density in 
Fig. 4(b). This is consistent with our analytical formula of the 
effective polaron model, Eq. (15), which does predict an in- 
duced phase below the critical phase 1/2 A^^ for the parameters 
used in Fig. 5. In Fig. 5(b) we have chosen the same parame- 
ters as in Fig. 5(a) but with a large temperature. The finite tem- 
perature in this plot manifests itself in a thermal background 
distribution superimposed on the coherent distribution, which 
reduces the contrast of the ToF image. However, even at such 
large temperatures the features of the distribution are still vis- 
ible and we expect this method to be able to distinguish be- 
tween the distributions with a central peak or a central dip at 
typical experimental temperatures. 

An application of this phase detection method via ToF im- 
ages might be to use the impurity as a nondestructive probe 
to reveal the rotational state of the underlying BEC. If the ro- 
tation is sufficiently close to a critical phase, the probe will 
acquire a phase which is detected by the central dip in the ToF 
image. 



IV. CLUSTERING OF IMPURITY ATOMS 

The presence of the coherent phonon cloud suiTounding the 
lattice atoms mediates a long-range interaction between dis- 
tant atoms. Furthermore, the onsite interaction decreases with 
the polaronic level shift Ep. For a sufficiently large mediated 
interaction this can result in an attractive onsite potential [24]. 
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FIG. 6: (Color online) (a) Correlations of impurities for Ng = 16 
lattice sites, 3 atoms, and Ua ~ —0.125.7a. The two lines represent 
the two families of correlations depending on the momentum state of 
the system, i.e., the phase a^. Each family is a set of correlations 
for many values of a^. (b) The spectrum exhibits crossings below 
and above the critical phases Ociit (the first critical phase is indicated 
with a vertical dashed line), where the impurities switch their mo- 
mentum ground state (thick solid and dashed lines). The thin gray 
lines indicate higher lying energy levels. 
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FIG. 7: (Color online) Cumulative level spacing for different disorder 
strengths e/ Jc = 0, 1, 20 (solid blue, dashed red, dashed-dotted yel- 
low lines, respectively). The edge of the dark shaded area indicates 
the Poisson prediction and the light shaded area the GOE prediction. 
The curves describe a crossover from regular behavior to chaotic, 
and back to regular for increasing disorder. This is for 3 BEC atoms 
and 2 impurities in lattices with 5 sites averaged over 100 indepen- 
dent disorder realizations. The other parameters are Uc = 1.5Jc, 

Ua = O.lJc, Ul = Jc, Uc = 0.057. 



phase term in the attractive Bose-Hubbard model is to intro- 
duce energy crossings and to bind the atoms more strongly for 
certain values of the phase. 



In this section, we will see how a phase twist influences the 
clustering properties of the impuiities. All computations in 
this section are based on numerically exact diagonalization of 
either the one- or two-species Bose-Hubbard Hamiltonian. 

We study the effective single-species Bose-Hubbard model, 
Eq. (4), which contains a phase term ~ aa in the hop- 
ping, and assume an attractive interaction iUaJ Ja < 0). The 
effect of the mediated interaction is measured by the density- 
density correlation {fijfiji) between sites j and j'. In Fig. 6(a) 
we see the effect of the phase term on the density-density cor- 
relations )■ Depending on the value of the induced 
phase the ground state is either in a weakly or a more strongly 
bound state. The strongly bound state with a higher onsite 
coiTelation is observed when the ground state of the system 
changes its momentum. In the energy spectrum this switch 
coiTesponds to a crossing of the two lowest energy levels. The 
energy spectrum of the system in Fig. 6(b) shows such cross- 
ings close to the first critical phase aait = 1/2 A^s. We find 
two such crossings below and above each critical phase. For 
the values of aa in between two such crossings the ground 
state becomes stronger bound. Thus the two families of cor- 
relations indicated in Fig. 6(a) consist of the coiTelations for 
all Qa either between two crossing above and below a criti- 
cal phase or outside the crossing region. For different sets of 
parameters the picture may feature more than two differently 
bound sets of states depending on the number of crossings 
within the lowest energy level. In conclusion, the effect of a 



V. QUANTUM CHAOS 

We now briefly address the question if our system exhibits 
quantum chaotic behavior as observed recently for a single- 
mode Bose-Hubbard model with a phase twist [29]. In our 
two-component Bose-Hubbard model, the phase of the second 
component is induced through coupling with the first com- 
ponent. This is different from the one-species model, where 
such an interspecies coupling is not present and where a phase 
term has to be assumed a priori. To formalize the notion of 
quantum chaos we denote with /S.Ej the distance between 
two neighboring energy levels Ej and Sj+i, and with IS.E 
the average over all spacings A£'j. Furthermore, we define a 
quasi-continuous parameter s = AE/AE for the normalized 
energy spacings. Then a nonchaotic system follows the Pois- 
son distribution with a probability density p{s) ~ cxp(— s). 
On the other hand, a quantum chaotic system follows a GOE 
distribution ]3(s) = 7rsexp(— 7rs'^/4)/2. The repulsion of lev- 
els expressed in the GOE is a result of correlations in the 
system, which are not present in the uncorrelated Poissonian 
(nonchaotic) statistics. 

Because of the various symmetries in the Bose-Hubbard 
model one cannot expect a global quantum chaotic behavior 
of such a system. One approach is to restrict oneself to the lo- 
cal behavior in suitable subspaces which do not exhibit these 
symmetries [27]. Another approach — and typically easier to 
achieve in experiments — is to break symmetries explicitly, 
for example, breaking translational symmetry by introducing 
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FIG. 8: (Color online) tests of the Poisson (solid line) and GOE 
predictions (dashed line). For intermediate values of u (0.2 < m < 
0.9) the energy distribution closely follows the GOE prediction. Be- 
low and above these values we observe a more regular behavior as 
the system is in either the hopping-only or the no-hopping state, re- 
spectively. As parameters we have chosen TVs = 5 lattice sites, 3 
BEG atoms, 2 impurities, Qc = 0.057. The results have been aver- 
aged over 100 independent disorder realizations. The inset shows the 
individual cumulative level distributions (thin lines, which appear as 
a gray area at this resolution) together with the Poisson (solid line) 
and GOE predictions (dashed line). 

random disorder [29]. We introduce disorder into the two- 
component Bose-Hubbard Hamiltonian with a phase twist by 
adding a local term with a random strength 

The energies Cj and are two independent random variables 
uniformly distributed on an interval (— e/2, e/2). In Fig. 7 we 
plot the integrated level distribution I{s) = p{s')ds' for 
different values of the disorder strength e. For vanishing dis- 
order, the system follows closely the regular behavior of the 
Poisson distribution. Choosing the same order of magnitude 
for all energy scales Ja,c — Ua,c — C// ~ e we recover a simi- 
lar behavior as observed in the one-component Bose-Hubbard 
model in [29]. The energy level distribution follows most 
closely the distribution of a GOE, which indicates quantum 
chaos. Increasing the disorder further reestablishes the regu- 
lar behavior because the atoms tend to localize in the disor- 
der potential. The effect of the interspecies interaction on the 
level spacing is investigated in Fig. 8. There we have defined 
a parameter u — Uj / e — Ua,c/e with Ja.c/e. — 1 — u. As 
we sweep u from to 1, the underlying Hamiltonian changes 
from describing an ideal two-component system with hop- 
ping to a fully interacting system without hopping. In the 
absence of the disorder potential these two extremal cases 
would lead to degenerate eigenenergies (going from momen- 
tum Fock eigenstates to spatial Fock states). Disorder lifts 
the degeneracy and we can observe global quantum chaotic 
behavior. We recover very similar curves with the definition 
u = Ui/e and Ja.cJ^ — Va,cl& = \ — u. The tests 
plotted in Fig. 8 indicate that the system initially follows a 
regular behavior but quickly changes to a GOE distribution. 



Increasing the interaction further, ultimately the behavior be- 
comes regular again in the no-hopping regime. We have de- 
fined = EJ[^num(sJ)-/(sJ)]^ where /nuni(sj) are the nu- 
merically obtained cumulative level spacings at point Sj and 
/(sj) is either the integrated Poisson or GOE distribution. 



VI. CONCLUSION 

We have presented a method for creating an artificial mag- 
netic field in a ring of trapped neutral atoms. Our method 
works by submerging atoms trapped in an optical lattice into 
a BEC which exhibits a phase term. We have discussed dif- 
ferent setups which lead to an induced phase in the impurities. 
First, with a BEC in a rotating ring, second, with a rotating 
BEC in a static ring. We have then derived an effective po- 
laron model, which reduces the problem to a single-species 
Hubbard model with a phase term on the hopping. For realis- 
tic parameters our analytical formulas predict induced phases 
up to aa — 0.03 in the first setup or \aa\ — 0.25 with a 
BEC in an excited rotational state. Comparisons with nu- 
merical solutions of the full problem for small systems show 
that our model can be used to qualitatively describe the sys- 
tem even beyond the constraints set by the analytical deriva- 
tion. Furthermore, we have discussed methods for observing 
the induced phase in the impurities: by using ToF images and 
measuring mass currents. The ToF images show a sharp tran- 
sition when a critical value of the induced phase is crossed. 
Increasing temperature or interspecies coupling obscures the 
sharp transition but the main features of the transition remain 
observable. Finally, we have compared and extended studies 
of quantum chaos to the two-species Bose-Hubbard Hamil- 
tonian with a phase twist. By introducing random disorder 
we could observe the onset of chaos when all energy scales 
in the system are of the same order, similar to the result of a 
one-species Bose-Hubbard Hamiltonian with a phase twist. In 
the two-species model the interspecies coupling ensures that 
the impurities exhibit a nontrivial phase, which is necessary to 
observe quantum chaos in this model. 

Contemporary optical lattice technology allows the inves- 
tigation of large systems of coherent neutral atoms evolving 
under a precisely known Hamiltonian, which is often not the 
case in actual condensed matter experiments. The methods 
presented here can be seen as a simulation of the effects of 
magnetic fields on ring systems. However, the ID treatment 
also allows us to gain insight into a 2D system of impurities 
immersed into a rotating 2D BEC [17]. Our results suggest 
that, in general, one may expect a decrease in the induced 
phase for nonvanishing interaction in 2D. However, the reduc- 
tion is not expected to destroy the large phases obtainable in 
2D. Furthermore, detection methods similar to the ones pre- 
sented here will be applicable in the 2D setup. In 2D a phase 
in the hopping term of the Bose-Hubbard Hamiltonian gives 
rise to the quantum Hall effect [5]. In light of this prospect, an 
extension of our methods to 2D will be of great interest and 
offer a quantum simulator for even more condensed matter 
models. 
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keeping terms up to second order in the fluctuations. For the 
second order part we keep two operators and replace the re- 
spective other two by Nq in the interaction terms of the Hamil- 
tonian. Furthermore, we substitute the chemical potential and 
use the bosonic commutation relation [d^, d^] = 1. This pro- 
cedure yields the Hamiltonian 



^(0+2) ^ 



It: 



90 



Ucna) 



Appendix: Bogoliubov approximation in tlie lattice with a phase 

twist 



n ( "-90+9 
"90-9 



To derive the Bogoliubov excitations for a BEC in a ID 
lattice with a phase term in the hopping we follow [54], where 
the authors derive the Bogoliubov approximation for a Bose- 
Hubbard model without a phase twist. First, we include 
the chemical potential /i^ in the Bose-Hubbard Hamiltonian, 
Eq.(ll), 



(A.l) 



We rewrite this Hamiltonian in the momentum represen- 



tation with 



Noting that 



^ gi27ri(9-g')/A'= ^ NsS^y we arrive at 



where 



D 



90+9 



UcUq 
Ucflo 



90-9 



Ucflo 



Here we have defined £qg±q ■— £q„±q ~ ^qg- The quadratic 
terms are diagonalized by the Bogoliubov transformation 



"90+9 

4o-9 



r> f'^qo+q 

Ul-9 



(A.2) 



In order for bq and 6^ to obey the bosonic commutation rela- 
tion the coefficients have to fulfill jugp — — 1. Plugging 
this expansion into the effective Hamiltonian and requiring 
that the resulting matrix {B^^Y DB^^ is diagonal we arrive 
at the following conditions for the coefficients 



H 



BH 



E 



A^c ) dldq 



_C4 
2N. 



E dldl,dg„dg,„Sq+q>^ 



',q"+q"' I 



q,q',q",q"' 



where iq = —2JcCos{27rq/Ns — 9) and 9 = 2TTac- For 
the case considered in Sec. 11 B the BEC is assumed to be in 
the ground state with integer winding number qo, which min- 
imizes the energy iq. We therefore write 9 = 2TTqo/Ns + 
A9/Ns, where A9 E [— vr, tt) determines the mismatch of the 
externally given phase twist 9 and the phase of the ground 
state. The ground state energy is then given by e^^ ~ 
-2JaCos{A9/Ns). 

The ground state is occupied with a macroscopically large 
number of atoms in mode qo, that is {(fq^^dq^^) ~ {dggd^qg) ^ 
and all other momentum state densities with g ^ are 
negUgible. This allows us to rewrite the creation and anni- 
hilation operators in terms of a mean field contribution and 
quantized fluctuations, that is 



d d^ 
dq,d\ 



d„,dl for q^qo. 



The chemical potential attains a value such that the Hamil- 
tonian is minimal with respect to Nq- It is given by /i^ = 
e^^ + UcfiQ, where no = Nq/Ns is the density of the con- 
densed atoms. An effective Hamiltonian is derived by only 



«~90+9l"9l 



^go-9l"9l +Ucno\Uq~Vq\ ^kuj^, 
'2eqUqVq + UcfloiUq - Vq)"^ = 0, 



where — (e 

9 ^ 

AJcSin^{TTq/N,)cos{A9/Ns). 
the solution 



90+9 



+ 

These 



90-9 



)/2 
equations 



have 



huji 



Ei 



1 = 



q 



- 1 



where we have defined Eq ~ ^ ^qi^q + 2Ucno) and := 

{li+q - ei-q)l2 = -2JcSm{27Tq/Ns)sm{A9/Ns). Thus 
the Hamiltonian up to second order is 



7j(o+2)=-ic/,f^g7V. + iE' K-(^-9 

9 

+ lT!^qHK- 



UcTlo)] 



For vanishing phase twist A^ = and the excitations have an 
energy given by Eq, which has the form of the free Bogoli- 
ubov energy Eq but with the free dispersion replaced by the 
dispersion relation of the first band in the lattice. For small q 
we see that Eq is linear in the quasimomentum, which means 
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(b) Here we have shifted the summation index by go in order for 
the Bogoliubov transformation, Eq. (A. 2), to be applicable. 
After application of the transformation we find 



H 



(0+1) 



j 



where the coupling between the impurities and the Bogoli- 
ubov modes is given by 



FIG. 9: (Color online) Contours of the relative condensate density 
no/n at r = for (a) iV^ = 10 and n = 0.5, and (b) Ns = 30 
and n = 2. The values in the plot indicate the relative density of the 
respective contour. 



that the quasiparticles behave like phonons as in the continu- 
ous case. For higher q the spectrum resembles massive parti- 
cles in a lattice. Turning on a phase twist results in an asym- 
metry in the spectrum because Kq is an odd function of q. 
Hence, the system of quasiparticles will prefer the branch of 
quasimomenta with lower energies, which will result in a drift 
of quasiparticles. 

Next we focus on the interaction Hamiltonian Hj between 
BEC and impurities, Eq. (12). As in the continuous case, 
we only keep terms linear in the fluctuations in the interac- 
tion Hamiltonian. Rewriting it in the momentum representa- 
tion and introducing a macroscopic population of one mode 
go yields 



3 

UiVNq 



'9+90 



j 9 
J2TTjq/N, 



-i2TTjq/Ns 




* \ - 



-aixjqjN^ 



The condensate density is defined as the macroscopic, 
condensed part of the total density n — N^/Ns. The noncon- 
densed part then arises from the expectation value of the fluc- 
tuation number operators, that is fi = uq + {d]^dq) /Ns- 

Substituting the fluctuations via Eq. (A. 2) we see that is 
implicitly given by 



' 9 



2E^ 



2 



1 



The Boltzmann factor enters through averaging over the 
phonon distribution at a temperature T. In Fig. 9 we plot the 
relative condensate density fiQ/n for two different filling fac- 
tors and lattice sizes. It becomes clear that our mean field 
approach is not valid close to the critical values of Uc at very 
low interaction Uc because the condensate density vanishes 
[see Fig. 9(a)]. On the other hand, we see that for small inter- 
action we can assume that ~ n provided that Uc is not too 
close to a critical value. 
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